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The forces on individual electrons in an unscreened nondegenerate electron fluid, due 

00 . 

to electron density fluctuations, have been calculated using Monte Carlo simulations 
and determined experimentally over a broad range of the plasma parameter T. The 
experimental results are obtained from the magnetoconductivity cr(B) measured for 
electrons on liquid helium below 1 K for B < 8 T. The magnitude and density 
dependence of <j{B) are explained by the many-electron theory of magnetotransport. 
The internal electric fields found from the experiments are in excellent agreement 

O ■ 

with the simulations. 
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Electrons above the surface of superfTuid helium often form a two-dimensional (2D) 



normal fluid M: for characteristic electron densities n ~ 10 12 m~ 2 the interelectron distance 
~ n~ 1 ' 2 exceeds the de Broglie wavelength, and the system is nondegenerate; at the same 
time, the ratio of the characteristic Coulomb energy of electron-electron interaction to the 
kinetic energy, the plasma parameter Y = e 2 ( y Tm) 1 ^ 2 / 'Ane^kT is large, and therefore there is 
short-range order in the electron system ||. For V > 127 (low T) electrons form a 2D crystal 
0. Mean-field effects such as long- wavelength plasma oscillations are well understood for 
a normal 2D electron liquid H but much less is known about the detailed behavior of 
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individual electrons. An experimental probe is required while, on the theoretical side, the 
problem is complicated by the absence of "good" quasiparticles. 

In this Letter, we present data from Monte Carlo simulations and experimental mea- 
surements of an important but unexplored characteristic of a 2D electron fluid, the internal 
electric field Ef that drives an electron as a result of its interaction with other electrons. 
Unlike the long- wavelength fluctuational fields known in plasma physics ||, the field Ef, 
although also of fluctuational origin, determines the force driving an individual particle, 
and is not described by the theory ||. The force on a particle is an important dynamical 
characteristic of a system. In "conventional" fluids with short-range interatomic inter- 
action, the forces have been a subject of extensive research Q. The mean square force 
remains finite in a fluid (though the mean square displacement diverges), but it would be 
expected to display singular behavior at the liquid-crystal transition. 

A special significance of the field Ef for a nondegenerate 2D electron system stems 
from the fact that, over a broad range of parameters, it strongly affects magnetotransport. 
In the single-particle approximation the electron energy spectrum in a magnetic field B 
perpendicular to the electron layer consists of discrete Landau levels, separation hu c (uj c = 
eB/m is the cyclotron frequency). The centers of the cyclotron orbits move only because 
of the random potential of the scatterers. Therefore scattering is always strong, and the 
pattern of transport is very different from the standard Drude picture which applies for 
weak coupling at B = and in which scattering events are short and well separated in 
time. 

The field Ef causes the centers of the cyclotron orbits to drift at a velocity Ef/B, and 
may "restore" the weak coupling picture. This effect was first discussed and investigated |7j 
for transport in the quantum limit Tiuj c / kT ^> 1 . It was shown later || Q that the field Ef 
also dramatically affects the magnetoconductivity cr(B) in classically strong magnetic fields, 
[iB 3> 1, huj c / kT < 1 (/i is the zero-field electron mobility), for vapor-atom scattering above 
1 K. In particular, it restores the Drude-type behavior of cr{B) for comparatively small B 
which has been known experimentally since fllOfl . Cyclotron resonance measurements |JTTJ] 
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have also demonstrated the importance of Coulomb interactions in this system. Several 
other mechanisms have been proposed to explain previous measurements of <j(B) at higher 



Our new measurements of er(-B) have been done below 1 K, in the ripplon scattering 
regime. In this range the mobility /i is extremely high (<2000 m 2 /Vs), enabling quantitative 
characterization of Ef as a function of density and temperature. 

We consider the distribution of Ef for a classical normal liquid; the results also apply 
for quantizing magnetic fields provided the motion of the centers of the cyclotron orbits is 
semiclassical. Since fluctuations in the system are thermal, and the field is due to electron- 
electron interaction, the scale for Ef is given by the characteristic field E : 

(E f 2 > = F(T)E 2 , E = (kTn 3 / 2 /4ne ) 1/2 . (1) 

The scaled dimensionless mean square field F(Y) can be easily found for large Y (low 
T) in the 2D crystal phase pi. Here, the force on an electron eEf arises because of 



the displacement of electrons from their lattice sites R n . In the harmonic approximation 
the force is linear in the displacements and has a Gaussian distribution. The function F 
incorporates contributions from both transverse and longitudinal modes of the crystal, and 
F(T) w 8.91, independent of T. 

In the opposite limit of a nearly ideal plasma, F < 1, the major contribution to the 
force on an electron comes from pair collisions, which gives F(T) w 2n 3 / 2 /Y. 

In the most interesting range of the electron liquid and the melting transition, the func- 
tion F(Y) was obtained from Monte Carlo simulations. We used the Metropolis algorithm 
and the Ewald summation technique as in 0(a), with periodic boundary conditions. The 
field on an electron was evaluated as the gradient of the potential in which the electron 
was moving. 

The results for Ef for the number of particles N = 196, 324 are very close to each other. 
In the range T > 30 the probability density distribution p of the field components E x , E y 
is close to Gaussian. The functions F(Y) and p are plotted in Fig. 1. The scaling function 
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F(T) decreases nearly monotonically with increasing V. However, remarkably, its variation 
is very small in the range T>10, although the structure of the system changes dramatically, 
from a crystal to a liquid with a correlation length of twice the mean electronic separation. 
The function F(T) has a smeared singularity at the melting point. Detailed discussion will 
be given elsewhere [111 . 



The magnetoconductivity of 2D electrons above superffuid helium was measured using 
4 mm diameter Corbino disk electrodes (see Fig. 2) 100 /xm beneath the electrons ||. Free 
electrons were held over the central drive electrode A, a ring electrode E, and receiving 
electrodes Bl, B2 and B3 surrounded by a planar guard G. An ac voltage Vq (typically 
10 mV) at a frequency up to 10 kHz was applied to electrode A and the ac currents / to 
the electrodes B measured. For a perfect conductor the phase of the capacitively coupled 
current / is 7r/2 with respect to Vo- The phase shift <f)(B) away from n/2 was measured 
for perpendicular magnetic fields B < 8 T for electron densities 0.5 x 10 12 <n<2 x 10 12 m -2 
at temperatures 0.6<T<0.9 K in the fluid phase. The phase shift 4>(B) is proportional to 
a^ 1 {B) for 0<O.3 rad, while for larger phase shifts the theoretical response function was 
used. The density n was determined from the -ve dc bias voltage on electrode E required 
to cut off the current between electrodes A and B. 

Figs. 2 and 3 show the measured magnetoconductivity for several densities and temper- 
atures. In low fields, B < 0.5 T, the data accurately follow the simple Drude-like result, 
even for values of fiB as large as 500: 

The electron mobility was determined from the B 2 dependence of cr~ l (B), as a function of 
density and temperature, and is in excellent agreement with previous B = measurements 



by Mehrotra et al. | 15| . The measured mobilities from 0.6 to 0.9 K are close to the 
theoretical values for a classical strongly correlated electron liquid |14[], with scattering 
by both ripplons and 4 He vapor atoms taken into account [pl|. The mobility is slightly 



density-dependent, primarily because the electric field that presses electrons against the 
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helium surface increases with n, and therefore so does the electron-ripplon coupling. The 
data is plotted as ne/fia(B) vs. B (Fig. 2) or B 2 (Fig. 3) using the experimental values 
of \i for each n and T. For B < 0.5 T the data lie on the universal line, ne/na(B) = B 2 
(line a in Figs. 2 and 3). 

At this point we should stress that the simple-minded Drude model (2) effectively 
applies because of the internal electric fields || [TJJ]. To understand the effect qualitatively, 
we notice that the many-electron system transfers its momentum to short-range scatterers 
via individual electron-scatterer collisions. In a certain range of the parameters, all that 
an electron "knows" about other electrons during a collision is the field Ef, and this field 
is time-independent if the collision is short enough. Then the Einstein relation for the 
conductivity applies: 



a(B) =ne 2 L 2 T-\B)/kT, (3) 

where L = L(B) is the diffusion length and r _1 (i?) is the relaxation rate. For fiB 3> 1 
the diffusion length is given by the mean radius R of the cyclotron orbit, L 2 = R 2 /2 = 
(h/2eB)(2n + 1), with n = l/[exp{hu c /kT) - 1]. 

It is the relaxation rate in (3) that is primarily affected by internal electric fields. 
For kT > eE { Xr> hu c (X T = h/{2mkT) 1 / 2 ) the field E f smears out the Landau levels 
and thus eliminates the effects of their discreteness. Therefore t^ 1 (B) = t _1 (0), and the 
^-dependence of a is given by that of L 2 = mkT/e 2 B 2 oc B~ 2 , i.e. the many-electron 
theory gives en/fia independent of n as observed, cf. the solid line a in Figs. 2 and 3. In 
contrast, in the single-electron picture the relaxation rate is increased by the density of 
states enhancement factor uj c t(B), since the states in the energy strip hu c are "compressed" 
down to the Landau level collision width hr^ 1 ^). The total magnetoconductivity from 
the self-consistent Born approximation (SCBA) for ripplon |16| and gas-atom scattering 
17| and their self-consistent combination cr s = (of s + ^gs) 1 '' 2 i s plotted (line b) in Fig. 2 



for n = 0.55 x 10 12 m~ 2 and T = 0.7 K. At 2 Tesla, the SCBA overestimates a(B) by an 
order of magnitude. 
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A distinctive feature of the classical cr(B) is saturation with increasing B. This arises 
because an electron in crossed Ef and B fields moves along a spiral with a step ~ 2nE{/ Buj c . 
The number of times it encounters a short range scatterer in the classically strong field 
B is then N cnc ~ Aj-i?co , c /27r£'f (Ay is the uncertainty in the position of an electron). For 
N enc > 1 one would expect the scattering rate, and thus cr(B), to increase by a factor 
~ A^nc oc B 2 , and we find (cf. Eq.(2)) 

ne „ „ f 2m 3 kT^ ' ~ 



mb)-'* Bl= V^) (E?>1/2 (4) 

The saturation of ne/fio~(B) is clearly seen in Figs. 2 and 3 for B > 1.0 T. The limiting 
value of ne/n<j(B) increases with density (Fig. 2) and temperature (Fig. 3) and is directly 
proportional to the rms internal electric field, from Eq.(4). It is this which enables the 
internal field to be determined experimentally. 

For B>1 T quantum effects become substantial. For hu c ^> kT the diffusion length 
L « (h/2eB) 1 / 2 , and also N enc ~ (H/eB) 1 / 2 Buj c /2irE i , so that ne/fia(B) oc B' 1 / 2 decreases 
with increasing B. For higher B>5 T (depending on n and T) the duration of a collision 
TcoI-ct(-B) and this theory no longer applies. 

Eq.(4) is written || for short-range scatterers (such as helium vapor atoms). In the 
case of scattering by ripplons, an extra factor arises which is numerically close to 1 for the 
experimental conditions here. Calculations of the magnetoconductivity for gas-atom and 
ripplon scattering were made using the full semiclassical many-electron theory in the range 
huu c >kT ||14|| , with values of (E 2 ) taken from Fig. 1. The results for the total many-electron 
magnetoconductivity a m = a Tm + a gm are shown in Fig. 2 (lines c, d and e, increasing 
n) and Fig. 3 (lines b, c, d and e, increasing T), and show satisfactory agreement with 
the experiments for t co \ t(B) (the extrapolation of the theoretical curves to the range 
r coi ~ T (B) is shown dashed). 

Conversely, the measured o~ l {B) at B = 2 T was used to obtain experimental values 
of the internal electric fields. This value of B is within the range of the applicability of 
the theory (kT > eE^h/eB) 1 ^ 2 ; htu c > kT) and is also far from the Drude region and the 
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region where collisional level broadening affects the scattering. The experimental values of 
(Ej?) 1 / 2 vs. E are shown in Fig. 4 (we renormalized E in Eq.(l) to allow for the dielectric 
constant of liquid helium). The points come from over 40 combinations of density and 
temperature between 0.6 and 0.9 K where the conductivity varies by more than an order 
of magnitude; no adjustable parameters have been used, and so the spread of the points 
may be considered to be within reasonable limits. Within the errors the measured field 
is oc Eq, with a constant of proportionality v — 3.11 ± 0.10. This can be compared with 
F 1 / 2 = 3.07 ± 0.03 from the Monte Carlo simulations for the range 20 < T < 70 covered 
by the experiments. A slight decrease of v with decreasing T seen in Fig. 4 lies within the 
errors. 

In conclusion, we have both computed and measured the internal electric fields in a 
nondegenerate 2D electron fluid, and the results are in excellent agreement. We show that, 
over a broad range of parameters, the magnitude and density dependence of the magne- 
to conductivity cr(B) of electrons on helium are determined by many-electron effects and 
can be understood qualitatively in terms of electron diffusion controlled by a fluctuational 
internal electric field. 
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Figure Captions 



Fig. 1. The scaled mean square field F(T) from Monte Carlo calculations. The asymptotic 
value of F for a harmonic Wigner crystal is shown dashed. Inset: the field component 
distribution. 

Fig. 2. ne/fia(B) versus B for n = 0.55 (o), 0.88 (y) and 1.89 (□) xl0 12 m- 2 at 0.7 K. The 
mobility // = 985, 830 and 520 ± 20m 2 /Vs, respectively. Solid lines a (low B) and c 
- e (high B) show many-electron calculations. Inset: Corbino electrode geometry. 

Fig. 3. ne/fia(B) versus B 2 for T = 0.6 (O), 0.7 (□), 0.8 (y) and 0.9 K(o). The mobility 
\x = 755, 620, 430 and 250 ±15m 2 /Vs, respectively. Lines a (low B) and b e (high 
B) show many-electron calculations. 

Fig. 4. Values of the internal field (E 2 ) 1 / 2 versus E at 0.6 (O), 0.7 (□), 0.8 (y) and 0.9 
K(o). Solid line shows the best linear fit with (E 2 ) 1 / 2 = 3.11E . 
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